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ABSTRACT 


The  purpose  of  this  contract  is  the  development  of  a  high-power 
broadband  amplifier  using  the  principles  of  a  traveling -wave  multiple- 
beam  klystron.  The  method  being  followed  is  to  (a)  design  the  beam  and 
the  interaction  impedances  by  a  computer  study,  (b)  build  a  single -beam 
prototype  for  checking  the  calculations  and  obtaining  final  data,  and 
(c)  construct  a  full-power  15-beam  tube. 

At  this  point  in  the  program,  the  beam  and  interaction  designs  have 
been  completed,  and  fabrication  of  the  single -beam  prototype  tube  is 
almost  complete.  Test  equipment  for  the  single-beam  tube  has  been  set 
up;  the  modulator  for  the  15 -beam  tube  hcis  been  designed,  and  major 
components  are  being  received.  The  electromagnet  for  the  15-beam  tube 
is  in  place,  and  the  single-beam  tube  is  to  be  tested  in  various  positions 
in  it. 


Conclusions:  All  basic  elements  for  the  single-beam  prototype  have 
been  designed,  fabricated,  and  tested.  The  assembly  is  almost  complete, 
and  processing  and  testing  of  this  tube  can  now  proceed.  No  major  prob¬ 
lems  are  foreseen. 
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REPORT 


A.  STATEMENT  OF  PROBLEM 

The  objective  of  this  program  is  to  design  and  build  a  C*lland  ampli¬ 
fier  of  high  peak  and  average  powers  and  broad  bandwidth,  using  sh* 
principle  of  a  multiple -beam  klystron.  The  objective  specification  is  in¬ 
cluded  as  Table  I. 

B.  BACKGROUND  AND  APPROACH  TO  THE  PROBLEM 

Following  earlier  efforts  to  build  traveling-wave  multiple-beam  kly¬ 
strons  in  which  the  beam  interactions  intermediate  between  the  input  and 
output  circuits  were  with  traveling  waves,  it  w^as  determined  that  the  re¬ 
quired  combination  of  bandwidth  and  efficiency  could  best  be  obtained  by 
the  use  of  traveling -wave  input  and  output  sections  w-ith  standing- w^ave 
(cavity)  intermediate  interactions. 

Based  on  these  findings,  the  approach  taken  in  the  present  work 
became: 

(1)  to  compute  a  suitable  set  of  interaction  impedances  for  a  single 
beam; 

(2)  to  compute  a  gun  design  to  supply  this  beam; 

(3)  to  build  beam  testers  as  necessary  to  demonstrate  the  beam 
design; 

(4)  to  build  a  single -beam  RF  tube  to  demonstrate  the  interaction 
design  and  obtain  optimum  impedance  values;  and 

(5)  to  build  a  multiple -beam  tube  incorporating  the  results  of  these 
preliminary  tests. 

Items  (1)  and  (2)  are  complete  except  for  refinement  as  a  result  of 
further  data.  This  report  covers  work  on  Items  (3),  (4),  and  (5),  arranged 
in  that  order. 
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no  time,  even  when  the  cathode  was  operated  temperature -limited,  was 
there  any  observable  change  in  pervearce  during  the  entire  100-micro- 
second  pulse.  Beam  transmission  through  to  the  collector  pole  piece  was 
better  than  99  percent,  with  a  magnetic  field  of  2550  gauss.  At  the  design 
magnetic  field  of  1830  gauss,  the  transmission  was  96.  7  percent,  with 
about  0.  9  percent  being  intercepted  at  the  collector  pole  piece  in  each 
case.  The  magnetic  field  profiles  corresponding  to  these  two  cases  are 
given  in  Figure  1.  Both  field  values  are  within  the  design  capability  of  the 
magnet  being  built  for  the  15-beam  tube. 

Figure  2  shows  the  cathode  current  as  a  function  of  voltage  for  vari¬ 
ous  values  of  cathode  temperature.  The  measured  perveance  of  1.27  is 
2  percent  below  the  design  value  of  1.  3  microamperes  volt^^.  The  cathode 
temperature  must  be  1030  C  to  avoid  temperature -limited  operation  at  60 
kilovolts. 

The  beam  tester  was  X-rayed  both  cold  and  with  the  cathode  at  full 
temperature  to  verify  the  estimates  of  thermal  expansion.  It  was  found 
that  the  expansion  of  the  entire  cathode  structure  was  0.  010  inch  greater 
than  had  been  estimated.  The  relative  position  of  the  cathode  with  respect 
to  the  focusing  cup  was  also  found  to  be  off,  but  only  by  0.  003  inch.  A  new 
beam  tester  will  be  constructed  with  the  correct  hot  dimensions  and  will 
be  tested  in  both  the  cylindrical  field  and  in  the  long  rectangular  field  being 
made  for  the  final  15-beam  tube.  Beam  performance  will  be  checked  in 
various  positions  in  the  rectangular  field  to  ensure  acceptable  operation  of 
the  multiple -beam  tube. 

Power  Dissipation  in  Loaded  Cavities 

Using  the  computer  data  for  design  No.  86  (imported  in  Quarterly 
Progress  Report  No.  3),  the  power  to  be  dissipated  in  each  loaded  cavity 
was  calculated,  and  the  results  are  plotted  in  Figure  3.  In  cavity  4  (the 
last  double-tuned  cavity)  the  peak  power  reaches  8000  watts  at  5500  MHz, 
due  to  the  sharp  rise  in  I/IQ  at  the  low  end  of  the  band  (see  Figure  5  in 
Quarterly  Progress  Report  No.  3),  while  in  cavity  5  (the  single-tuned 
loaded  cavity)  the  peak  value  is  8500  watts  at  5850  MHz,  due  in  this  case 
to  the  peaking  of  the  cavity  impedance  (see  Figure  4  in  Quarterly  Progress 
Report  No.  3). 

At  the  maximum  duty  of  0.  01,  these  two  cavities  must  therefore  be 
capable  of  absorbing  80  to  100  watts  average  power  without  becoming 
overheated. 
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PEAK  CURRENT -AMPERES 
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The  heat  conduction  equation  was  solved  for  a  rectangular  block, 
cooled  on  the  rear  face,  with  energy  entering  by  the  front  face  and  assumed 
to  be  absorbed  in  the  material  at  a  rate  which  decreases  linearly  to  zero 
at  the  back  face;  this  leads  to  the  following  equation  for  the  steady- state 
temperature  rise  at  the  front  surface: 


T 


2aW 

3AK 


(1) 


where:  A  is  the  cross-sectional  area  perpendicular  to  the  heat 

flow  (CM2) 

a  is  the  depth  parallel  to  the  heat  flow  (CM) 

K  is  the  thermal  conductivity  (watts/CM^/°C/CM) 

W  is  the  power  to  be  absorbed  (watts),  and 
T  is  the  temperature  rise  (°C). 

The  greater  precision  of  an  exponential  assumption  about  the  absorp¬ 
tion  is  not  justified,  because  of  considerable  uncertainty  in  the  value  of  K. 

Using  the  book  value  of  k\  0.  04  watt/CM^/°C/CM,  for  the  planned 
absorber  arrangement  in  cavity  4,  the  calculated  temperature  rise  was 
960°C.  While  this  was  obviously  excessive,  it  was  still  thought  to  be  an 
underestimate,  because  the  book  value  of  K  was  probably  applicable  to 
compact  SiC  of  near  theoretical  density,  rather  than  the  porous  material 
being  used. 

To  obtain  directly  relevant  data,  the  specific  heat  and  thermal  con¬ 
ductivity  of  an  actual  specimen  were  measured.  The  values  obtained  were: 

Specific  heat:  0.25  calories/GM 

Thermal  conductivity  K:  0.025  watt/CM^/  C/CM 

This  indicated  that  the  average  temperature  rise  would  be  even  worse 
than  first  calculated,  and  that  the  absorber  must  be  redesigned  to  decrease 
'a'  or  increase  'A',  or  both,  while  still  appearing  as  a  good  RF  match. 


1.  Handbook  of  Thermophysical  properties  of  Solid  Materials,  Goldsmith 
ct-al.  Vol.  3,  VU-D. 
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bp  (])  iS  perhaps  deceptively  simple.  With  a  fixed  power  W  to 

«  iJ*.  f?  >e  fiC,meS8  'a'  raiSeS  the  PO--  density  i/the 

material  at  the  same  time  it  shortens  the  conduction  path  to  the  cooled  wall 
Intuition  might  suggest  that  the  two  effects  would  cancel  out,  but  this  is  not  ’ 

the  a  f  impr°Ved  conduction  outweighs  the  higher  power  density,  so  that 
the  surface  temperature  becomes  proportional  to  the  thickness  One  should 
therefore  use  the  thinnest  layer  that  is  still  6  should 

RF  match.  7  11  capable  of  Providing  a  good 

the  U  AJher  t^alS’  tHe  arrangement  of  Figure  4  was  adopted/  Since 

sprina  Ht  coT  V  Plate  ^  inStalled  in  the  cavity  by  means  of  a 

spring  fit,  cooling  is  obtained  at  the  ends  of  both  arms  as  well  as  at  the 

center  where  the  piece  is  f  tened.  The  thicknegd  25  CM  Ud  the 

effective  area  'A'  is  9  CM*.  Inserting  these  values  in  equation  H)  the 

toP30C°C  S5Lr^  rir  18  appr°ximately  60°C;  if  the  cavity  wall  is  held 
to  JU  C,  and  the  SiC  surface  should  reach  approximately  90°C. 

average  ^  abs°rber  lay°ut  of  Fig^e  4  was  then  tested  at  full 

g  p  wer  by  exciting  it  with  a  C-Band  voltage  tunable  magnetron  and 
a  matched  coupling  loop,  inserted  through  the  beam  tunnel.  Stripes  of 

idTar:/rrenSi,iVe  We"  appUed  *°  «“  •”*«•  of  the  slucon  car- 

bide,  and  the  power  was  raised  to  80  watts.  Color  changes  in  the  paint 
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Figure  4  -  Final  Design  of  Loading  for  Cavities  2,  3,  and  4 
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indicated  temperatures  between  95  and  145°C  on  the  side  walls,  and  under 
95  C  on  the  thicker  blocks  at  the  back.  After  operation  at  120  watts  --  50 
percent  above  the  maximum  in-service  power  level  --  the  145°C  color 
change  occurred  over  the  central  area  of  each  side  piece,  indicating  a 
maximum  temperature  of  probably  160  to  170°C.  At  this  temperature,  no 
outgassing  problems  are  anticipated,  after  a  600°C  bakeout. 

The  experiment  indicated  that  this  design  is  adequate  for  the  average 
power  that  is  expected;  it  also  confirmed  that  the  assumptions  on  which 
equation  (1)  is  based  are  good  enough  for  this  kind  of  calculation,  in  which 
a  precise  result  is  not  needed  --  only  an  assurance  that  the  working  tem¬ 
perature  will  be  well  under  the  bakeout  temperature. 

Applying  the  same  methods  to  the  cavity  5  design,  the  calculated 
working  temperature  was  under  200  C,  and  no  change  in  this  design  was 
necessary. 


The  pulse  temperature  rise  was  also  calculated;  owing  to  the  deep 
penetration  of  the  power,  and  the  relatively  high  specific  heat  (2-1/2  times 
that  of  copper),  the  pulse  component  was  found  to  be  completely  negligible 
compared  to  the  average  temperature  rise,  for  both  cavities. 

The  dissipations  in  cavities  2  and  3  are,  of  course,  small  compared 
to  that  in  cavity  4,  so  these  pose  no  problem.  The  same  loading  is  used, 
however,  for  the  sake  of  uniformity. 

Cavity  and  Waveguide  Response  Curves 

When  the  revised  cavity  loading  had  been  tested  to  the  full  average 
power,  as  described  in  the  previous  section,  permanent  loads  were  fitted 
to  cavities  2,  3,  4,  and  5.  The  response  curves  of  these  cavities  are  shown 
in  Figure  5.  The  differences  between  these  curves  and  those  shown  in 
Figure  4  of  Quarterly  Progress  Report  No.  3  are  mainly  due  to  the  fact 
that  the  computed  curves  include  the  effect  of  beam  loading,  while  the  cold- 
test  experimental  curves  do  not.  The  responses  obtained  conform  closely 
to  the  desired  cold  frequencies  and  Q  values. 

The  response  of  cavity  6  is  not  shown,  since  without  beam  loading  or 
internal  loading,  it  is  simply  a  spike  at  6100  MHz.  The  frequency  of  the 
cavity,  as  brazed,  was  within  2  MHz  of  the  design  value. 
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Figure  5  -  Cavity  Response  Curves  and  Waveguide  Reflectometer 
Curves  (Traced  from  Oscilloscope  Photos) 
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Also  shown  in  Figure  5  arc  the  reflectometer  traces  of  the  two 
waveguides,  after  brazing  and  final  adjustment.  The  bandpass  is  shifted 
about  20  MHz  downwards,  compared  to  the  brass  cold-test  models, 
probably  due  to  small  differences  in  the  shapes  of  the  fixed  copper  tunnel 
noses,  compared  to  the  earlier  adjustable  ones.  The  VSWR  is  a  little 
high  at  the  high  frequency  end  of  the  band  as  a  result,  but  is  acceptable 
for  the  single-beam  tube.  The  multiple-beam  waveguide  design  is  being 
adjusted  to  recenter  the  bandpass  correctly. 

Single-Beam  Tube  Assembly 

Assembly  of  the  single-beam  tube  has  proceeded  more  slowly  than 
expected.  The  major  problem  here  has  been  distortion  of  the  waveguides 
adjacent  to  the  windows  during  brazing.  Several  windows  were  lost  during 
at  empts  to  rebraze;  the  alumina-ceramvar  braze  is  adequately  matched 
for  many  heat  cycles  up  to  600  C,  but  rebrazing  at  1000°C  takes  it  well 
above  the  kink  in  the  ceramvar  expansion  curve,  and  cracks  have  developed. 

ht  basic  cause  of  the  distortion  appears  to  be  the  use  of  standard  OFHC 
copper  waveguide,  which  was  perhaps  a  false  economy.  It  is  planned  to 
redesign  the  windows  for  the  15-beam  tube  so  that  the  ceramic  braze  is 
only  subjected  to  one  brazing  cycle  --  the  next  joint  being  a  weld  --  and  to 
machine  the  waveguide  sections  from  solid  material. 

At  this  time,  the  input  and  output  waveguides  and  all  the  cavities  have 
been  brazed,  aligned,  and  welded  together.  The  cathode  parts  are  all 
available,  but  are  not  being  assembled  to  the  tube  until  the  window  problems 
are  solved,  in  order  to  protect  the  cathode  surface.  Four  windows  are  on 
hand,  but  two  of  these  are  marginal,  and  completion  of  the  tube  is  being 
e  ayed  a  few  days  while  efforts  are  made  to  complete  two  fully  satisfactory 

winHnwo  ' 


Figure  6  depicts  the  tube,  with  the  body  section  welded,  the  external 
supporting  structure  in  place,  and  the  windows  and  cathode  insulator 
clamped  on. 

Test  Equipment 

The  test  position  for  operating  the  single-beam  tube,  and  providing 
the  physical  structure  and  magnetic  field  for  the  15-beam  tube,  was  com¬ 
pleted.  Measurements  on  the  magnetic  field  are  described  in  the  next 
section  (see  "Magnet"). 
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Figure  6  -  Single-Beam  Tube  Showing  Body  Section,  External 
Supporting  Structure,  Windows,  and  Cathode 
Insulator 


The  driver  tube  was  received  from  IT&T,  after  a  demonstration  run 
in  their  plant.  The  gain  of  this  tube  is  well  over  the  specified  minimum, 
with  the  result  that  a  100-mW  output  sweeping  oscillator  can  be  used  as  the 
primary  RF  source;  the  calibrated  frequency  sweep  and  the  levelled  output 
will  be  major  convenience  factors.  The  modulator  for  the  driver  tube  has 
been  installed  and  tes^d. 

The  20  by  20  foot  cage  for  the  15-MW,  150-KW  modulator  was  com¬ 
pleted;  all  major  components  are  on  order,  and  some  have  been  delivered. 

Magnot 


The  elongated  magnet  coils  were  received  late  in  the  reporting  period. 
The  winding  irregularities  were  a  little  greater  than  expected,  so  that  for 
the  single-beam  tube  it  has  only  been  possible  to  stack  8  instead  of  the  9 
coils  that  were  planned;  for  the  15-beam  tube,  which  will  have  no  tuners, 
there  is  no  restriction  requiring  that  the  gaps  be  lined  up  with  the  tuning 
screws,  and  for  this  tube  it  will  be  possible  to  use  9  (or  even  10)  coils. 
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However,  the  measured  magnetic  field,  as  a  function  of  is@ffcl  g#$l 
current,  is  within  1  percent  of  the  calculated  value.  Hence,  the  margin# 
of  safety  that  were  allowed  In  designing  the  coil*  are  net  needed,  and  it 
will  be  possible  to  operate  to  over  2500  gauta  with  8  coil*  while  still  r* « 
maining  within  the  design  maximum  current  for  each  ceil.  The  intended 
range  of  operation  is  1830  to  2500  gauss. 

With  plain  pole-plates,  the  measured  field  dips  6.  5  percent  at  the 
center  of  the  magnet  gap.  The  addition  of  1/2  by  1/4- inch  field  sir *ifh(«o«f * 
along  the  outer  edges  of  the  pole-plates  reduced  the  field  variation  acros# 
the  gap  to  less  than  2  percent,  as  shown  In  Figure  ?. 


End  effects  begin  to  be  measurable  who  i  the  beam  hole*  are  within 
3.  5  inches  of  the  end  of  the  pole-plate  structure.  Since  the  waveguide 
matching  sections  extend  at  least  3  Inches  from  the  end  beams,  the  pole* 
piece  extensions  will  contribute  at  most  one  inch  tc  the  length  of  the  tube 
body. 

Impedance  Taper  for  Output  Waveguide 

The  output  filter  circuit  of  a  hybrid  multiple-beam  klystron  consists 
of  two  regions  --  a  uniform-impedance  build-up  portion  in  which  the  voltage 
at  each  successive  beam  increases,  and  a  power  region  in  which  the  gap 
voltage  at  each  beam  is  held  essentially  constant  at  its  optimum  value  by 
reducing  the  circuit  impedance  at  each  beam  to  compensate  for  the  additional 
power  (see  pages  1-3  of  Quarterly  Progress  Report  No.  1).  Since  the  RF 
current  induced  in  the  output  cavity  of  a  klystron  is  a  function  of  the  RF 
voltage,  each  build-up  beam  will  induce  a  different  RF  current  into  the  out¬ 
put  eircuit.  Figure  8  shows  the  computed  relationship  between  current  and 
voltage,  plus  a  linear  approximation  which  is  suitable  for  voltages  below 
the  desired  power-  beam  voltage. 

Figure  9(a)  is  a  schematic  diagram  of  a  build-up  section  having  a 
uniform  impedance  ZQ.  Only  the  forward-traveling  current  waves,  If,  will 
be  considered  in  the  design  since  backward  waves  induced  by  the  beams 
will  tend  to  cancel  one  another.  Thus,  since  the  forward-wave  current  is 
in-phase  with  the  induced  current,  the  forward  currents  entering  and  leaving 
the  section  will  be  related  by: 

Jf(n+1)  =  Jf(n)  +  i  Xn+1  <2> 
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Figure  8  -  Normalized  RF  Beam  Current  Versus  Normalized  RF  Beam  Voltage 
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The  induced  current,  I^j,  may  be  related  to  the  gap  voltage  by  the  follow¬ 
ing  linear  approximation: 


I 


n+1 


(3) 


where  Ip  and  Vp  are  the  desired  current  and  voltage  of  the  power  beams 
and  Ip  £  the  zero  voltage  intercept  of  the  linear  approximation  as  shown 
in  Figure  8.  Solving  these  equations  for  the  voltage  at  the  (n+  1)  build-up 
beam  yields: 
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These  recursion  relations  allow  the  voltage  and  current  at  each 
build-up  beam  to  be  calculated.  When  the  equations  yield  a  voltage  exceed¬ 
ing  V  ,  that  beam  is  in  the  power  region  instead  of  the  build-up  region, 
and  it*s  impedance  must  be  reduced.  Figure  9(b)  depicts  a  power  section 
in  which  the  ideal  transformer  provides  a  matched  transition  between  the 
impedances  of  adjacent  sections.  In  this  case: 


^(n+l)  Jf(n) 


V  -  V  =  T  7 

p  n+1  f(n+l)  n+1 


Thus, 
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This  quadratic  equation  can  be  solved  to  give  Zn+i  in  terms  of  Zp  and 
^(n+l)’ 
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A  BASIC  computer  program  (THAL21)  for  calculating  the  mid-band 
characteristic  impedance  for  each  section  is  given  in  Appendix  A.  A  con¬ 
servative  preliminary  taper  design,  based  on  ninety  percent  of  the  calcul¬ 
ated  current,  is  given  in  Table  III.  A  final  design  can  be  computed  after 
the  current-voltage  response  of  the  single -beam  prototype  is  measured. 


Table  III  -  Preliminary  Taper  Design 


IQ  =  17.46  Amps  V  =  53, 550  Volts 

P 

I  =  14. 73  Amps  Z  =  2, 260  Ohms 

r  O 


Beam  Position 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Circuit  Impedance 
2260 
2260 
2260 
1470 
1069 
837 
686 
581 
503 
444 
397 
359 
327 
301 
278 


Gap  Voltage  (KV) 
18.66 
36.  30 

52.  98 

53.  55 
53.  55 
53.  55 
53.  55 
53.  55 
53.  55 
53.  55 
53.  55 
53.  55 
53.  55 
53.  55 
53.  55 
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Low-Impedance  Circuit 

For  the  taper  design  given  in  Table  III,  the  impedance  must  be  re- 
uced  by  a  factor  of  8.  1.  If  the  capacitively  loaded  waveguide  ured  in  the 
ui  d-up  section  were  scaled,  the  waveguide  impedance  would  have  to  be 
reduced  by  a  factor  of  8.  1,  and  the  gap  capacitance,  increased  in  the 
same  ratio.  The  resulting  waveguide  would  be  0.  0b2 -inch  high,  with  an 
interaction  gap  of  the  order  of  0.  020  inch.  These  close  spacings  would 
c  eate  mechanical  problems  and  prohibitively  high  RF  electric -field 
intensities,  leading  to  field  emission  and  arcing. 

An  alternate  configuration  for  the  low-impedance  sections  has  been 
investigated  It  consists  simply  of  a  uniform-width  waveguide  with  two 
regmns  of  different  heights,  as  shown  in  Figure  10(a).  This  circuit 
matches  the  phase  shift  and  impedance  response  of  the  capacitively  loaded 

.  .vegmde  very  ciosely,  but  has  a  gap  (H2,)*  of  0.  080  inch  for  the  lowest 
impedance  section. 

The  two-height  filter  can  be  accurately  designed  (neglecting  the 
beam  tunnel)  without  any  experimental  measurements.  Since  there  are  no 
ariations  in  the  transverse  direction,  the  zero-mode  cutoff  of  the  filter 
occurs  when  the  width  is  exactly  one-half  of  a  free-space  wavelength,  in 

fLld  pTtteTnatVb  a8/  Unif°rm-height  octangular  waveguide.  The  electric 
field  pattern  at  this  frequency,  wo,  is  that  obtained  electrostatically  bv 

applying  a  voltage  between  the  top  and  bottom  faces  (with  the  side  walls 
removed). 

Th!  71  CUt°ff  °r  resonance  short-circuit  symmetry  planes 

the  ends  of  the  section  shown  in  Figure  10(a).  This  i  -mode  cavity  is 
eqmvaient  to  a  portion  of  ridged- waveguide  having  a  length  equal  to  the 

of  th1  W’d  \  W  ’-a  M  ^  7r"m°de  resonance,  w  ,  W1  must  equal  one-half 
of  the  ridged -guide  wavelength.  Thus,  the  ridged-guide  cutoff  frequency, 
0>r  is  given  by:  H  y' 

/  T~<*>  z' 

W  =  0)  l/1  -  ,a. 

r  7T  V  &  (o) 


Variables  occurring  in  the  computer  programs  are  written  thus,  rather 

than  in  subscript  notation  (H2),  to  facilitate  comparison  of  the  equation,, 
with  the  programs.  ^ 
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Figure  10  -  Physical  Layout  and  Equivalent  Circuits  for 
Output  Waveguide  Sections 
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The  ridged-guide  cutoff  may  be  calculated  from  the  resonance  of  two 
lenoths  of  TEM  line  (having  widths  of  W  1/2  plus  a  discontinuity  capacitance, 
Cd  ,  as  indicated  in  Figure  10(b).  That  is: 


CO  C  , 
r  d 


W1 

754H1 
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W1 

754H2 


tan 


CO  L2] 
r 

c 


0 


(9) 


where  c  is  the  velocity  of  light.  The  filter  design  is  simplified  by  assum 

ing  the  following  approximate  relation  for  C,: 
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1.  1 


(picofarads/meter)  for 


Hi<5 

H2 


(1C) 


and 


If  C0Q  and  u)jj.  (from  which  cor  and  W1  are  determined)  plus  LI,  L2, 
one  of  the  heights  are  specified,  the  approximation  allows  direct  solution 
for  the  other  height.  If  HI  is  specified,  the  resulting  equation  for  H2  is 
linear,  whereas  specifying  H2  yields  a  quadratic  in  HI. 


The  phase  shift  and  impedance  response  of  the  filter  may  be  calculated 
from  the  equivalent  circuit  shown  in  Figure  10(c).  At  n  mode,  this  equiva¬ 
lent  circuit  yields: 


W 

d  '  [754H1J 


+ 


tan 


(11) 


Because  of  the  relationship  between  the  three  frequencies  (to  ,  co  ,  and 
Ct>r),  the  arguments  of  the  tangent  and  cotangent  functions  in  ftiis  equation 
and  the  one  given  previously  for  the  ridged-guide  are  the  same.  Multiply¬ 
ing  the  ridged-guide  equation  by  the  radical  and  subtracting  yields: 
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2.  J.  R.  Whinnery  and  H.  W.  Jamieson,  "Equivalent  Circuits  for  Dis¬ 
continuities  in  Transmission  Lines",  Proc.  IRE,  Vol.  32,  pp  98-114 
Feb.  1944. 
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Since  LI  and  L2  may  be  adjusted  to  tune  to  any  frequency  without  alter¬ 
ing  the  value  of  C^,  in  the  previous  equation  may  be  taken  as  any  arbi¬ 
trary  frequency,  CO.  Thus,  the  discontinuity  susceptance  is  a  parallel 
resonant  circuit  with  a  resonant  frequency  equal  to  the  guide  cutoff  fre¬ 
quency  and  a  capacitance  equal  to  the  discontinuity  capacitance  of  TEM 
lines  with  the  same  heights  and  one-half  the  width. 


The  phase  shift  per  section  may  be  calculated  from  the  open  and  short- 
circuited  admittances  of  a  half-section  as  indicated  in  Figure  11  (also  refer 
to  page  6  of  the  Quarterly  Progress  Report  No.  1). 


/3L 


2  tan 


2  tan 
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The  characteristic  impedance  at  the  plane  of  the  beam  is: 


Zq  (beam) 


1 

Vy  Y 
oc  sc 


(14) 


The  characteristic  impedance  at  the  plane  midway  between  beams  is: 


Z  (mid)  = 
o 


(15) 


A  BASIC  program  (THAL22)  for  designing  and  analyzing  these  filters 
is  given  in  Appendix  B.  The  required  input  quantities  are  fG,  f n  ,  HI  or 
H2,  LI,  and  L2.  The  program  computes  the  other  height  and  the  width  plus 
0  L,  Z0(beam),  and  ZQ(mid)  across  a  specified  frequency  range. 

Figure  12  shows  the  computed  phase  shift  response  for  a  filter  of 
this  type  and  the  experimental  resonances  of  a  four-section  model  with 
and  without  the  beam  holes.  A  capacitively  loaded  waveguide  designed  for 
the  same  beam-to-beam  separation  and  cutoff  frequencies  has  a  computed 
response  that  agrees  within  a  phase  shift  of  about  one  degree  at  all 
frequencies. 


In  order  to  maintain  good  coupling  to  the  electron  beam,  the  height 
(H2)  must  not  exceed  0.  170  inch.  Furthermore,  L2  must  exceed  0.  100  inch 
to  allow  for  the  0.  100-inch- radius  beam  hole.  As  a  result,  the  two-height 
circuit  is  most  suitable  for  impedance  levels  below  700  ohms. 
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Figure  12  -  Frequency-Phase  Shift  («-/3  )  Plot  for  Output  Waveguide  (Power  Section) 


Intermediate -Impedance  Circuit 


Intermediate  values  of  interaction  impedance  (i.e.  between  700  and 
2200  ohms)  can  be  obtained  with  a  hybrid  circuit  employing  the  two  differ¬ 
ent  height  regions  defined  in  Figure  10  plus  capacitive  tunnel  noses.  Since 
this  circuit  has  one  more  degree  of  freedom  than  the  two-height  filter,  it 
is  possible  to  specify  both  heights  (HI  and  H2)  as  well  as  both  lengths  (LI 
and  L2)  and  the  two  cutoff  frequencies.  The  dependent  variables  are  the 
width  and  the  gap  capacitance,  Cg.  A  direct  solution  is  not  possible,  but 
an  iterative  computer  solution  is  convenient. 

At  77  mode,  a  voltage  minimupi  (short  circuit)  exists  between  beams, 
and  a  maximum  (open  circuit)  at  the  beams.  Thus,  referring  to  Figure  11 
and  adding  Cg/2  at  the  left  side  of  each  circuit,  we  obtain: 
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The  zero  mode  has  open  circuits  at  both  planes  so  that: 
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The  admittances,  Yac  and  YQC,  are  not  simply  related  to  the  width,  Wl, 
or  the  waveguide  cutoff  frequency,  £q.  However,  the  problem  may  be 
linearized  by  assuming  that  for  small  variations  in  fc: 
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Thus,  for  an  initial  estimate  of  fc,  we  have: 
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respectively,  the  required  width  is  0.983  inch.  The  resulting  phase  shift 
and  impedances  at  the  center  and  edges  of  the  operating  band  are: 


Freq  (MH,) 

1 3  L  (degrees) 

ZQ  (beam) 

ZQ  (mid) 

5370 

43.  8 

3802. 

1430. 

5650 

102.6  ^ 

2260. 

535. 

5030 

152.8 

3899. 

192. 

Design  of  Tapered  Output  Circuit 

A  preliminary  design  for  a  tapered  output  circuit  has  been  accom¬ 
plished  using  the  circuit  forms  and  computer  programs  described.  Table 
IV  gives  the  dimensions  and  mid-band  performance  of  each  section.  The 
design  of  the  high-impedance  build-up  circuit  had  already  been  determined 
by  general  considerations  such  as  cathode  diameter,  beam  diameter,  and 
beam  coupling  as  described  in  earlier  reports.  Although  a  certain  amount 
of  computer  cut-and-try  was  involved,  the  following  general  steps  were 
followed  in  the  taper  design: 

1.  Determine  the  value  of  ZQ  (beam)  required  for  each  section 
from  the  RF  voltage -cur rent  characteristics  of  the  beam 
using  the  program  of  Appendix  A. 

2.  Investigate  geometries  of  the  two-height  circuit  that  are  suit¬ 
able  for  the  lowest  impedance  (last)  section  using  the  BASIC 
program  of  Appendix  B.  The  impedance  is  essentially  directly 
proportional  to  H2  and  decreases  as  L2  increases.  Check  that 
HI  is  reasonable  and  that  the  RF  electric  field  intensity  at  the 
band  edges  (where  the  impedance  is  highest)  is  not  excessive. 

3.  Investigate  geometries  in  which  H2  equals  0.  170  inch  (the 
maximum  gap  for  good  electron  coupling)  and  HI  equals  0.  500 
inch  (the  height  of  the  uniform  circuit),  to  determine  the 
maxim. '.m  impedance  level  and  thus  the  lowest  beam  position 
index  for  which  the  two-height  circuit  is  suitable.  In  this 
case,  section  8  is  the  first  of  the  two-height  design. 

4.  Construct  a  smooth  taper  of  values  of  HI  for  the  two -height 
sections  since  abrupt  changes  in  HI  represent  an  unwanted 
discontinuity  between  sections.  (Note  vhe  progression  of  HI 

in  section  8  through  15  of  Table  IV.)  For  each  of  these  sections, 
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Equivalent  length"  of  0.  3309  used  in  calculation* 


adjust  LI  (keeping  Li  plus  L2  constant)  and  interpolate  to 
obtain  the  desired  impedance  level  (if  possible).  Note  that 
the  ratio  of  ZQ  (beam)  to  ZQ  (mid)  at  mid-band  stays  constant 
(at  1.65  in  this  caje)  withip  a  few  percent  so  that  there  is  a 
cold  mismatch  between  each  section  equal  to  the  relative 
change  in  ZD(beam). 

5.  Construct  a  table  of  values  of  ZQ  (mid)  for  the  hybrid  stages 
to  provide  a  smooth  transition  between  the  build-up  sections 
and  the  first  two-height  section  (e.g.  let  ZQ  (mid)  vary 
exponentially  in  the  hybrid  sections).  (Note  the  progression 
of  ZQ  (mid)  in  sections  4  through  7  of  Table  IV.  )  Since  the 
hybrid  sections  have  one  more  degree  of  freedom  than  the 
two-height  ones,  the  ratio  between  ZQ  (beam)  and  ZQ  (mid) 
may  be  varied. 

6.  For  each  hybrid  stage,  let  HI  equal  0.  500  inch  (the  height  of 
the  uniform  circuit)  and  vary  H2  and  LI  to  obtain  the  required 
values  of  Zq  (beam)  and  ZQ(mid).  Note  that,  in  general  terms, 
changing  LI  changes  the  value  of  both  impedances,  whereas 
varying  H2  alters  the  ratio. 

The  dimensions  given  in  Ta.ble  IV  have  been  chosen  to  yield  identical 
computed  cutoff  frequencies  for  all  sections.  However,  there  is  a  slight 
discrepancy  between  the  calculated  and  computed  frequencies,  as  illustrated 
in  Figures  12  and  13.  These  deviations  result  from  neglecting  the  beam 
hole  in  the  two-height  circuit  and  the  finite  dimensions  of  the  capacitive  gap 
in  the  hybrid  circuit.  Cold-test  models  of  a  few  different  sections  will 
yield  estimates  of  the  shift  in  cutoff  frequencies,  and  the  sections  can  be 
redesigned  accordingly.  For  example,  Figure  12  shows  that  the  beam 
holes  increase  the  zero-  and  ir  -mode  frequencies  23  and  37  MHz,  respec¬ 
tively.  Thus,  the  sections  which  are  similar  to  this  model  should  be  re¬ 
designed  by  specifying  cutoff  frequencies  too  low  by  the  respective  amounts. 
In  the  case  of  section  number  eight,  this  redesign  changes  H2  from  0.  157 
to  0.  159  inch,  and  the  width  from  1.  114  to  1.  118  inch,  with  the  other 
dimensions  unchanged. 

The  two-height  circuit  lends  itself  to  an  alternate  impedance  tapering 
scheme  which  may  be  advantageous  where  a  larger  transformation  ratio  or 
a  higher  power  capability  is  required.  In  this  technique  the  beam  hole  is 
moved  transversely  from  the  center  of  the  guide  towards  the  side  wall 
where  the  electric  field  is  weaker.  Since  the  electric  field  distribution  is 
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sinusoidal,  the  impedance  transformation,  [L  , 
by: 


fJ- 
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2 

COB 


(~  n  D 
wi 


ia  given  (for  small  holes) 


(23) 


where  D  is  the  displacement  from  guide  center.  In  this  type  of  circuit  the 
hole  represents  a  relatively  minor  perturbation  so  that  moving  it  towards 
the  side  should  have  only  a  small  effect  on  the  phase  characteristic. 


D.  CONCLUSIONS 


1.  The  major  components  of  the  tube  --  cathode,  gun  structure, 
input  and  output  waveguides,  loaded  cavities  and  collector  --  have  been 
built  and  tested  as  close  to  full-power  operating  conditions  as  possible. 

2.  Assembly  of  these  components  into  a  single -beam  RF  tube  is 
almost  complete. 

3.  Test  equipment  for  the  single -beam  tube  is  complete. 

4.  The  magnet  for  the  15-beam  tube  is  complete. 

5.  The  modulator  for  the  15-beam  tube  is  in  the  early  stager  of 
construction. 

6.  The  two  most  difficult  theoretical  parts  of  the  program  --  the 
beam  interaction  calculation  and  the  output  waveguide  impedance  taper  -- 
are  complete. 

7.  No  insurmountable  technical  problems  are  foreseen  in  the  attain¬ 
ment  of  the  objectives  of  this  project. 

E.  RECOMMENDATIONS  AND  FUTURE  PLANS 

1.  The  immediate  objective  is  the  completion  and  testing  of  the 
single-beam  tube,. 

2.  Parts  for  the  15-beam  tube  which  do  not  depend  critically  on 
the  single-beam  data  will  be  machined. 


3.  The  beam  interaction  and  impedance  taper  calculations  will  be 
re-run  using  measured  values  from  the  single-beam  tests,  and  drawings 
for  the  cavities  and  output  waveguide  will  be  revised, 

4.  When  testing  of  the  single-beam  tube  is  complete,  the  test 
engineers  will  work  on  construction  of  the  150-KW  modulator. 
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Appendix  A  -  Computer  Program  (THAL21)  for  Calculating 
Impedance  Taper 


100  DATA  19«4# I 6» 37# 53550# 2260 
110  REN  IMPEDANCE  TAPER  CALCULATION 
120  READ  I0#1I#VI#Z0 
I  30  REM  AMPS#  VOLTS#  UNO  OHMS 

1  40  PRINT  M INDUCED  I"#"INOUCED  2"#'*V  MAX"#"Z<0>" 

150  PRINT  "AT  VO"  #"AT  V*V  MAX"#"  "#  "UNIFORM" 

160  PRINT  I0#I!#V1#Z0 

170  PRINT  "BEAM  NO."#"Z<N>"#"Z<N>/ZCO>"#"I  FWD"#"V<N>" 

180  LET  SaSOR(ZO) 

190  LET  1"0 

200  REM  CALC  VOLTAGE  <V>  AND  FORWARD  CURRENT  CI>  IN  BUILD-UP 
BIO  FOR  Mai  TO  15  STEP  1 

220  LET  VZO*<  1*10/2  >/<  1  -•  5*Z0*C I  I-I0>/V1> 

230  REM  IF  VOLTAGE  EXCEEDS  OPTIMUM  CALCULATE  REDUCED  Z 

2  40  IF  V»VI  THEN  290 

250  LET  I-l«(I0*<II-I0>*V/Vl>/2 
260  PRINT  M.ZO."  1"#I#V 
2  70  NEXT  M 

280  REM  CALC  I  AND  LINE  IMPEDANCE  AT  EACH  POWER  BEAM 

290  FOR  N-M  TO  IS  STEP  I 

300  LET  Aall/2 

310  LET  BaI*S 

320  LET  Ca-Vl 

330  L£T  R9«S0R(B*B-4*A*C> 

340  LET  S1»<-B*R9>/<2*A> 

350  LET  Z>S1*S1 
360  LET  IaI*S/SI*l 1/8 
370  PRINT  N#Z#Z/ZO# I 
380  LET  S-SI 
390  NEXT  N 

400  LET  Pa. 5E-6*Z*I*1 
410  LET  B"2*VI/I l/ZO 

420  PRINT  "PWR  OUT  <MW)"#"EOUIV  Z/.5Z<0>" 

430  PRINT  P#B 
4  40  PRINT 
450  GOTO  120 
460  END 


34 


Appendix  B  -  Computer  Program  (THAL22)  for  Designing 
Output  Waveguide 


100  OAT A  5300*5990.0.. OB.. I5».  55 

101  DATA  5300. 5**0. .3180***0..  IS*.  55 

110  ACM  DESIGN  OF  TV0-HEI6HT  WAVEGUIDE  FILTER 
I  80  RCAO  rO.F9.HI.H2. LI* L2 
I  30  REN  MHZ  AND  INCHES 

140  RfcH  HI  IS  UNKNOWN  UN c CSS  H2«0  IN  WHICH  CASE  H«  IS  IMMMOIM 
I  SO  LET  WI«5*0I • 4/FO 

ISO  FRINT  MF(0>"."F(P1  )"."LCNBTN(  1  >"."LCNBTH< •>"» " VI DTH~ 

I  TO  FRINT  F0.F9.LI.L2.WI 

ISO  REN  CALCULATE  REQUIRED  F<CUT0FF)«F3  OF  TRANSVERSE  R1DGE0  GUIDE 
1*0  LET  F3«r**S0R<l-<ro/r*>»2> 

200  REM  CALCULATE  TAN(BL)  FOR  LI  ANO  L2  AT  P3 
210  LET  W-6.263IBES*F3 
>20  LET  V>l « I802*EI0 
230  LET  T|aTAN(W*Ll/V) 

240  LET  T2«TAN(WaL2/ V) 

250  REM  SOLVE  EQUATION  FOR  HEIGHT  HI  (QUADRATIC)  OR  Ht  (LINEAR) 

2 SO  LET  Aa 10. 5334E* I2*W 
>70  IF  H2a0  THIN  330 
200  LET  Ba-H2«T2/A 
2*0  LET  C“*H2/TI/A 
300  LET  R*a5QR(B*B-4*C) 

310  LET  Hla(.B*R*)/2 
320  00  TO  340 

330  LET  H2a(HlaHI*A«T2*HI)/(HI*A«!/TI> 

340  LET  Z|aTS4*HI/WI 
350  LET  Z2a754*H2/Wl 
360  LET  FI-FO 

370  LET  C-B2. 445*(H I/H2- I )/FI 

3 BO  PRINT  "HEIOHT(l)"#"H2IGHT(2>"»"C(DlS>"»"F(CUTOrr>" 

390  PRINT  HI.H2.C.FI 

400  PRINT  "FREQ". "BETA  L"»"Z(0)  AT  BEAN"."Z(0)  MID" 

410  FOR  F-S370  TO  5*30  STEF  140 
420  REM  CALCULATE  WAVEGUIOt  FARAMETERS 
430  LET  R0a|-(PI/F)f2 
440  LET  R|aSQR(ABS(RO)) 

450  LET  B«5.3235E-44FaRI 
460  REM  RAOIANS  FER  INCH 
470  LET  Sa.| 

480  LET  T|aTAN(B*LI) 

490  LET  T2-TAN(B*L2) 

500  HEM  CALCULATE  X*S  OF  LINE  I  OPEN  (XI >  AND  SHORTED  (Yl> 

510  LET  X|aZiaT|/RI 
520  LET  Y|aS*ZI/TI/RI 

530  REM  AOO  DISCONTINUITY  SUSCEFTANCE  (B9> 

540  LET  B9a6< 2B3 1 BE-6*F*C 
550  LET  B9«B9*R0 ' 

560  LET  X|aXI/(l-XI*B9> 

570  LET  YlaYI/(l-YI*B9> 

5 BO  REM  TRANSFORM  X»S  ALONO  LINE  2 

590  LET  Z3aZ2/RI 

600  LET  X4aZ3*(XI*Z3*T2>/(Z3t3*XI*T2) 

610  LET  Y4aZ3*(Yl*Z3*T2)/(Z3*S*T l*T2) 

S20  REM  CHECK  WHETHER  FILTER  PROPAGATES 
630  IF  (*X4/Y4)  »0  THEN  STO 
640  PRINT  F 
6  SO  GOTO  T90 

660  REM  CALCULATE  BETA  LAND  Z(0> *S  AT  BEAN  AND  MID  PLANES 
670  LET  02* 57. 29SB*ATN( SQR ( *X4P  Y4>  > 

6 BO  LET  D2a2*02 
690  LET  ZaSQR(*X4*Y4) 

700  LET  Z4aZI/RI 
710  LET  X2aZ3*T2 
720  LET  Y2aZ3/(S*T2) 

730  LET  X2aX2/ ( I •X2SB9 ) 

740  LET  Y2aY2/(l«Y2*B9> 

750  LET  XSaZ4*(X2*Z4*TI)/(Z4«S*X2*TI) 

760  LET  YSaZ4*(Y2«Z4*Ti)/(Z4«S*Y2*TI) 

770  LET  YaSQR(>XS*YS) 

780  PRINT  F.02.Z.Y 
790  NEXT  F 
BOO  PRINT 
BIO  GO  TO  120 
820  END 
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Appendix  C  -  Computer  Program  (THAL23)  for  Designing 
Capacitively  Loaded  Waveguide 


1 00  DATA  5300  •  9990*  •  S»  •  34*  •  38*  •  38 

I  10  REM  C-L0A0E0#  TV0-3TEP  WAVEGUIDE  OESIOM  AMO  CALCULATION 
I  80  RCAO  F0«F9*HI#H2*LI«L8 
130  REM  MHZ  Mr  iNCHES 
I  40  PRINT  "F<ZEPJ>"»"F(P1>" 

I  50  PRINT  ro#r» 

160  PRINT  "HEIGHT!  I  >"*"HEI0HT<2)"*"LENGTH<  I  >"*"LEN0TH(2>" 

ITO  PRINT  HI«H8<  LI*L8 

ISO  LET  P8a(P04P?)/2 

190  REM  START  OP  ITERATIVE  LOOP 

800  POR  rial  TO  7  STEP  I 

810  REM  CALCULATE.  B< TEND) aRl 

880  LET  PaPO 

830  LET  PI >pe 

8  40  00 SUB  550 

850  LET  BBa.Riaviarus 

860  GOSUB  750 

870  LET  BI«B8 

280  REM  INCREMENT  P(CUTOPP)  ANO  RECALCULATE  B(ZER0)aB3 

290  LET  P|a|.02aP2 

300  GOSUB  550 

310  LET  BBa-RiaviaTlaS 

320  GOSUB  750 

330  LET  BSaBB 

340  REM  CALCULATE  B(PI»aB2 

350  LET  PaP9 

360  LET  Plare 

370  GOSUB  550 

380  LET  B8a-H|«Y|/Tl 

390  GOSUB  750 

400  LET  B2-BB 

410  REM  INCREMENT  P<CUTOPP»  ANO  RECALCULATE  B(PI >aB4 

480  LET  P|a|.084P8 

430  GOSUB  550 

440  LET  B8a-R|aYl/TI 

450  GOSUB  750 

460  LET  B4aB8 

470  REM  CALCULATE  DC'DFCC)  AT  ZERO  ANO  PI  MOOES 
480  LET  00* ( B3-B II/!. 0L*F2 > 

490  LET i 09a(B4-B8)/ (•02*P2) 

500  REM  CALCULATE  NEW  ESTIMATE  OP  P(CUTOPP> 

510  LET  P8aP8«<F94BI-P04Be>/(P0ab9«P94D0> 

520  NEXT  N 
530  GOTO  780 

540  REM  SUBROUTINE  TO  CALCULATE  HAVE6U10E  PARAMETERS 
550  LET  R0-I-(P1/P;r8 
560  LET  RlaSGR(ABS(RO>) 

570  LET  Ba 5. 323SE* 4*PaR| 

580  REM  RADIANS/ INCH 

590  LET  Ca88. 445*(HI/H8-1 l/P I 

600  LET  HlaS901.4/Pl 

610  LET  Y|aWl/TS4/HI 

680  LET  Y2aWi/7S4/H8 

630  IP  F*FI  THEN  700 

640  LET  Sal 

650  LET  ElaEXP<e*8*LI> 

660  LET  E2aEXP(8*B*L8> 

670  LET  Tla(El«l)/(EI61> 

680  LET  T2a(E2>l )/(82*l ) 

690  GOTO  730 
700  LET  Sa.l 
710  LET  TlaTAN(8*LI) 

780  LET  T8aTAN(B*L2) 

730  RETURN 
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740  REM  SUBROUTINE  TO  ADD  DISC  8  AND  TRANSFORM  B  AL0N8  LINE  S 

790  LET  B0aB8*6«283l6£”6*F*C*RO 

760  LET  8BaRI*Y24(88-S*Y2*Nl»TB>/(Y2*RI-BB*T2> 

770  RETURN 

780  LET  Cl— BI/<64283I8E-6*P0> 

790  LET  Wla990l.4/F2 

800  PRINT  "C<aAP>/B"«,•WIDTH••*MrCCUTOrr>•• 

810  PRINT  CI*WI«P2>.I06*(H2-.I7) 

820  LET  ZI»794*HI/WI 
830  LET  Z2a7S4*H2/Wl 
840  LET  P|a990l.4/Wl 

890  PRINT  "PRESt,«MGCTA  L"*"Z(0>  AT  BEAMM*MZ<0>  MID" 

C  60  FOR  F*9690  TO  9691  STEP  100 
870  Q0SU8  990 

680  REM  CALCULATE  X'S  OF  LINE  I  OPEN  (XI)  AND  SHORTED  (Yl) 

890  LET  X |aZI 4TI/RI 
900  LET  YUSaZ|/TI/RI 

910  REM  ADD  DISCONTINUITY  SUSCEPTJWCE  (D9> 

920  LET  B8a6»2S3l8E*6*F*Cl 
930  LET  89a6.283l8E-6*F*C 
940  LET  B9a89«R0 
990  LET  X|aX|/(|-XlaB9> 

960  LET  YlaYI/(l-YlaB9) 

970  REM  TRANSFORM  X*S  ALONG  LINE  2  AND  ADD  QAP  SUSCEPT4MCE  (M> 

o an  |  it  71a70/RI 

990  LET  X4aZ3«(XI*Z34T2>/(Z3»S*X|4Ta> 

1000  LET  X4aX4/(|oX4*aa> 

1010  LET  Y4aZ3a(YI«Z34T2>/(Z3*S*VlaTB> 

1020  LET  Y4-YV(l-V4aB8> 

1030  REM  CHECK  WHETHER  FILTER  PROPAGATES 
1040  IF  <*X4/Y4>  »0  THEN  1080 
1090  PRINT  P 
1060  OOTO  1200 

1070  REM  CALCULATE  BETA  L  AND  ZCOI’S  AT  BEAM  AND  MID  PLANES 
1080  LET  0Ba97.2998aATN(SGR(-X4/Y4>> 

1090  LET  DBaBaDB 
I  100  LET  ZaS0R(-X4*Y4) 

I  1 10  LET  Z4aZI/RI 
I  120  LET  X2'Z3*T2 

1190  LET  YBaC34<>l«B8*Z3*TB>/(Z34BB*8«Tt) 

I  140  LET  X|aXC/(l*X2*B9) 

1190  LET  V2aYf£S<|.YB*09> 

1160  LET  X9aX4*(X8*Z4aTI >/<(4*f*A2»TI > 

1170  LET  Y9«Z  *•(  Y2*I4*TI  )/<Z6*S*YI*TI) 

I  ISO  LET  YaS0ll(»X9*V9) 

I  190  PRINT  r»te»Z«Y«E/Y 
1200  NEXT  P 
1210  PRINT 
1220  OOTO  120 
1230  (MO 
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